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ABSTRACT. We report a quantitative comparison between theoretical predic- 
tions and observations of the intensity of galactic cosmic rays near the 
interplanetary current sheet. Comparison of our model calculations is 
made with a statistical analysis of observations of galactic cosmic rays 
at Earth and the simultaneous position of the current sheet. We use an 
ensemble of different current sheet inclinations, in order to make the 
analysis of the computations approximate the method used to analyse the 
data. 

1. INTRODUCTION. The transport of cosmic rays in the interplanetary 
magnetic field i6 the consequence of four basic effects - diffusion, 
convection, cooling and gradient and curvature drifts. The resulting 
transport equation for the distribution function f may be written: 

£ ^ v-(g- -y.pf v-v £ r u) 

where P is momentum, £ is wind velocity and ft* is the diffusion tensor. 

A variety of solutions of this equation^Tias appeared in the litera- 
ture. A straightforward application of the equation leads to a situation 
in which the particle drifts play a very important and perhaps dominant 
role in the modulation of galactic cosmic rays (see, e.g. Jokipii and 
Kopriva, 1979, Jokipii and Davila, 1981, and Kota and Jokipii, 1982). 
The problem of the intensity near the interplanetary current sheet was 
first addressed observationally by Newkirk and Lockwood (1981). Subse- 
quently, Kota and Jokipii (1982) showed that the Newkirk and Lockwood 
results were consistent with their three-dimensional code. 

More-sophisticated analyses of the observations were reported by 
Newkirk and Fisk (1985) and Newkirk, Asbridge and Lockwood (1985). Here 
we determine the agreement of these observations with our model. 

2. The Model Calculations . We have developed a 3-dimensional code to 
solve the full transport equation. We use a straightforward extrapola- 
tion of observations taken near the ecliptic plane. The model is static 
in a coordinate frame corotating with the Sun. The numerical technique 
and details of the model are described elsewhere (Kota and Jokipii 1982). 
The interplanetary magnetic field is the same as that used by Jokipii and 
Thomas (1981). The field at the solar surface is assumed to be uniform 
and radial, with opposite signs on either side of a magnetic equator. 
This magnetic equator is a tilted plane at the sun, which is a wavy 
neutral sheet at larger heliocentric radii in interplanetary space (see 
Figure 2 of Jokipii and Thomas 1981). The degree of waviness increases 
from a minimum near solar minimum to a large value near 90 near sunspot 
maximum. The field on either side of the current sheet is an Archimedean 
spiral. The case in which the northern hemisphere field is outward 
corresponds to the field parameter A being positive. The solar wind 
speed was taken radial and constant in magnitude at 400 km s~*. We are 
not yet able to use a spatially-varying solar wind speed. 

The parallel diffusion coefficient, Kii » was assumed to be inversely 
proportional to the magnetic field strength, B , 

K, - ( B eartt/ B ) <» 

with F being the particle rigidity in GV, is the particle velocity in 
units of velocity of light, and K 0 is a normalization constant in the 
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range 10^1 - 10^3 if is expressed in cm^/sec (see, eg*, Jokipii and 
Davila, 1981). The ratio of perpendicular and parallel diffusion coeffi- 
cients was kept constant at K Ct/^n “ 0*05-0.10. To reduce computing 
requirements the outer boundary was set at r = 15 AU, with some runs 
carried out for a 30 A*U. boundary. The general features of the 
solutions have been discussed in detail elsewhere (Kota and Jokipii, 
1982), and their discussion will not be repeated here, except to re- 
iterate that the solutions are clearly affected by the particle drifts. 
Of most interest here is the dependence of the intensity on the structure 
of the interplanetary current sheet illustrated in the contour plot in 
figure (1). It is clear that the magnetic field organizes the cosmic-ray 
intensity in a characteristic manner relative to the current sheet, and 
that this organization is different for the two signs of the interplane- 
tary magnetic field. Note also, for future reference, that the intensity 
is not simply dependent on the distance from the current sheet. 


Fig. 1 Computed contours of 
equal intensity at a radius 
of 1 A.U., for protons of 
energy 2.36 GeV. Inclina- 
tion of current sheet ■ 45 °. 
KqHL. 0 x 10 ^ cm^/sec, and the 
outer boundary was at 15 
A.U. All other parameters 
as in Jokipii and Kota, 
1982. 
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3. Comparison With Observations 

a. The Data. 

Given the above, it is of interest to study observat ionally the 
relationship between the current sheet and cosmic rays. The analyses 
carried out by Newkirk and Lockwood (1981) and Newkirk and Fisk (1985a) 
and Newkirk, Asbridge and Lockwood (1985) attacked this problem statisti- 
cally, in order to minimize the problem of transient, time-dependent 
effects. They studied the dependence of the intensity of the Mt. 
Washington neutron monitor (to protons of about 5.3 GeV energy) on the 
distance from the current sheet, for data obtained in the years 1973- 
1978. Data from the years around sunspot maximum were not included. The 
determined the position of the current sheet from coronal white light 
data, which has been determined to be quite accurate. Their figure (9) 
shows a scatter plot of the daily intensity vs. the heliomagnetic lati- 
tude at the point of observation (defined as the angular distance in 
degrees along a meridian to the current sheet, which defines the helioma- 
gnetic equator). The data show considerable scatter, which is to be 
expected since effects of transient disturbances, etc. may be expected to 
disrupt the pattern. Nontheless, there is a clear trend in the data, and 
the solid line gives the best fit of the function 

I ■ a Q + sin 2 ( ^ 5 ) (3) 

to the data, where ^*^is heliomagnetic latitude. The fit to the data, 
with an =2Ho7 and aj = - U m 7 is quite well-determined. In a subsequent 
paper, TJewkirk, Asbridge and Lockwood examined the energy dependence of 

the ratio aWag by analysing data from a variety of other sources. These 
data show that the effect depends inversely on energy, scaling approxi- 
mately as T~ 0 * 8 . It should be noted, however, that this obser- 

vational result is still consistent with there being no instantaneous 
latitudinal gradient. For, since all of the observations are taken near 
the solar equatorial plane, the large values of occurred when the 
current sheet was far from the equatorial plane, and the inclination of 
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»| 1 «K C i U i rrent sbeeb wa ® large. Hence, if the intensity of cosmic rays 
gl obally were small when the current-sheet inclination was large, then a 

i s 3 observed ®?iT 8 n^ 6 ,?- 8 *^ wo H ld tend ?° be lower at large m2 , as 

it doe. n« affec ” che thlS '»«>■•* h-l, .. 

b* The Model Calculations. 

In order to compare these data with the model, it is important to 
simulate the methods used to analyse the data very accurately. 

•*5 8t, k i We note that the inclination of the current sheet varied 
nft? S K derably °^® r th ® time Period spanned by the data set, so it would 
• e appropriate to use only one simulation, with one current sheet 
vflfiin at ^° n ’ tb ® tbe Efl rth carries the point of obser- 

t or S 1 ly 86Ven degre r es abov ? and below the heliographic equa- 

of m 1 • ^ th .® “ a * ln ! um value of mg m the data set shown in -figure 4 (9) 
nftion of 4?o nd Fl ?J “ approximately 55°, a maximum current sheet incli- 

Thomas flORH «JXi2j lve a P a . xlmun >. of 52°. The work of Jokipii and 

Thomas (1981) suggests a minimum lnclfiiation of the order of 15° Hence 

ree different inclinations, 15, 30, and 45 degrees were computed for 
f. C . h p P«T ter S6t * • i , f ? urbh inclination of 0 degrees was also used in 
* ua ?» * ltb . 1 f tt . le change in the results. Then, to simulate 
nnmW^f 5 he E £* fc u ln lt8 orbit * the intensity was determined at a 
P ° la u tS ^ wh i cb ware . spaced in heliographic latitude and longi- 
tude just as the Earth is in its orbit. For each such point, in addition 

SUfeknXfifkTaf cil&i-atS! »s‘ s »y 
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Fig. 2 a. Scatter plot 
of intensity vs /L m& 
summed over three incli® 
nations 15, 30 and 45 

degrees, for the para- 
meters of fig 1. b,c,d 
give the individual con- 
tributions for the three 
inclinations. 


^ (DE6) ^ 

r V U Vui n ® 8 ?* °^. P°.* nts ^ or three inclinations and 

along the Earth's orbit is illustrated in figure (2a) for one set of 
parameters. Note that there is considerable scatter in the points, 
reflecting the fact that the computed intensity depends on other para- 
5 s we ^ 1 . as A mg* . Nontheless, the calculated values show a trend 

obtained for the individual inclinations. 

Finally, a least squares fit of equation (3) to the synthesized 
scatter plot was made to obtain the "theoretical value" of the coeffi- 
a k for each given set of parameters. This procedure was 
f® p ® abed {°\ e . aca set of parameters (diffusion coefficient normalization, 
energy, etc). As mentioned above, the absolute values are irrelevant, so 

fu! l a ^?l' a k i8 u* 8 ® d i 1 ? w £ at foll ows. Illustrated in figure (3) gives 
izatTonfl Vnr rat *° for a range of diffusion-coefficient nolmal- 

lzations for particles of energy 2.36 GeV, together with the value for 
this energy given by Newkirk, Asbridge and Lockwood (1985). Clearly, 
there is a broad range of plausible diffusion coefficients for which the 
computed value is close to that observed. 
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Fig. 3 Dependence of computed 

a t /ao on K n for the parameters in 
figures 1 and 2. Horizontal band 
is observed value. 


We also studied the energy dependence of the ratio a^/ao*. We fo “ 
that although the calculated value does indeed decrease with increasing 
particle energy, and agrees well with the data at a few GeV energy, the 
functional dependence is not the simple power law seen in the data. We 
expect that changing the energy dependence of the diffusion tensor couJ. 
improve agreement here, but have not been able to verify this. 

4. Summary and Interpretation The analysis presented above demonstrates 
that a model of modulation, in which the dominant physical effect is the 
large-scale structure of the interplanetary magnetic field, is quantita- 
tively consistent with the analysis of the data reported by Newkirk and 
Lockwood (1981) and Newkirk and Fisk (1985a, b). At present the only 
significant discrepancy appears to be in the somewhat different energy 
dependences obtained from the data and from the model, aud it appears 
likely that modifying the diffusion coefficient may improve this. Taken 
together with other comparisons of the theory with data, including t 
prediction of 22-year solar magnetic cycle effects (Jokipii and Thomas, 
1981, and Kota and Jokipii, 1982), and detailed comparison with the 
inclination of the current sheet (Smith and Thomas, 1985), this suggests 
that many features of the solar-cycle modulation of galactic cosmic rays 
are a consequence of the model, and that drifts may well be a dominant 
process in solar modulation. 
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